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MMC #538 - OCEAN WATER COLOR ASSESSMENT 
FROM ERTS-1 RBV & MSS IMAGERY 
D. S. ROSS PR175 

PROGRESS REPORT, PERIOD 8 MAY - 7 JUL 1973 
CONTRACT NAS5-21862 

1. ■ Objective . The objective of the contract is to determine 

whether or not it is feasible, by Interactive image processing, 
to separate image densities representing ocean scene reflec- 
tances in the 460-490 nm region, from the full 460-580 nm 
band recorded by ERTS-1 RBV-1. By so doing, an indication 
of blue and green water color ratios may be obtained as an 
indication of fertility for fish populations. The method 
employed is based on photographically masking the RBV-1 
image with an MSS-4 image of opposite polarity. As the 
MSS-4 records a 490-600 nm band, subsidiary masking with 
RBV-2 ( 580-680 nm) is required to remove image densities 
attributable to 58 O- 6 OO nm scene reflectances. Additionally, 
masking to reduce residual radiometric correction errors is 
used. 

2. Summary of Work During Period . Radiometric error correction 
masking on the RBV-1 and -2 Images was completed, and these 
Images were remade at higher gammas to offset atmospheric 


2 



scene contrast degradation effects. A new image was made 
by combining RBV-2 pos and MSS-4 neg at matching gamma, to 
cancel image densities in the 580-600 nm spectral region. 

The resulting image appears to show this effect in the 
shallow water areas, but density reproduction control for 
the next step was Judged to be Inadequate. Re-calculation 
of atmospheric haze effects and the less than y 1.0 sensor 
reproduction of the original in the dark areas shows that 
scene reproduction gammas in the Images to be worked with 
require additional corrections. Work has been re-started 
and different methods of gamma and density control are being 
employed for Improved control. 

3, Conformance with Work Schedule . Delays in receiving Imagery, 
beyond the control of GSPC and this investigator entailed the 
submission of a new Data Handling Plan, which was approved 
17 Apr 73. Current work is on schedule, but the delay in 
starting will cause completion of the project to be in 
December 1973, rather than August 1973. A request to extend 
completion date, at no additional cost to the Government, 
has been submitted to the Technical Monitor. 
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Progress 


4, Correction masks for reducing residual radiometric shading 
errors were made as described previously, and applied to 
RBV-1 anJ -2 images of ERTS-1 image set 1007-15165, 30 Jul 12 

I 

(Caribbean area) for generating corrected negatives. Estima- 
tion of scene contrast reduction, due to atmospheric effects, 
was made by measuring image densities of light and dark 
subjects (Cloud and water areas) for gamma correction in 
reproducing the RBV-1 and MSS-4 Images to contrast ranges 
more nearly representing their true ground level albedo. 
Initially, RBV-1 and -2 were reproduced at y 3.85 as 
positives, and MSS-4 as a negative at the same y. 

5, RBV-2 pos and MSS-4 neg were combined in registration and 

a masked positive was made at y 1.0, to preserve the contrast 
of the y 3.85 working images. The objective was to produce 
a new image containing water image densities in the 490-580 nm 
range. The MSS-4 image covers a 490-600 nm spectral scene 
range, and the RBV-2 a 58 O- 68 O nm range. The clear, calm 
ocean water surface makes little or no spectral contribution 
to the image above 580 nm, and the water surface in this 
spectral region would be near Dmax in the RBV-2 positive, 
and clear in the negative; when combined with the 490-600 nm 
MSS-4 positive, scene reflectances above 580 nm would be 
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masked on the land areas, or where light in very shallow 
water is reflected from the bottom; but the surface of the 
deep water would print through. The masking effect was 
found in a shallow area between North and South Bimini 
Islands, and is shown in Figure 1. 

/' ■ 

Problems j 

6, Great difficulty . was experienced in maintaining exposure 

and density control to ensure the proper level of masking, 
although use was made of cloud and water image areas for 
reference, and related densities in the GSFC 15 step 
tablet. The steps in the latter and their density Increments 
do not lend themselves to convenient sensitometric control. 
Images of the same clouds, deep water and shallow areas in 
each spectral band, were used for exposure control in 
reproducing the range of image densities concerned on the 
straight line part of the reproduction film, in each case. 
Because high gammas were entailed, straight line image 
densities became very high when negatives and positives 
were combined for reproduction, net densities being above 
5,0 in some cases, and control based on image densities 
alone eventually became undependable. 
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Figure 1 (Photographic Print) 

ERTS 1007-15165, 30 Jul 1972 

Left . Print made by masking RBV -2 positive image 
with MSS -4 negative at controlled gammas and 
densities. A 20 nm yellow spectral region of 
bottom reflectance in the water has been suppressed 
between North and South Bimini Islands, in the lower 
part of picture. The same effect is found in the 
shallow water north of Grand Bahama Island, at the 
top of the picture. 

Right , Unmasked MSS-i| image for comparison. 
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It also became apparent that the image densities of 
primary importance, the ocean surface, were being affected 
during reproduction by their departure from y 1.0 in the 

I 

original j and work with this set of reproductions was 
halted to reexamine the problem. 

8. The RBV images match each other for registration, but do not 
match the MSS-4 image. A best fit is made in the north-south 
central parts in registering the RBV and MSS Images, since this 
is the main area of deep and shallow water under examination. 
Images corrected for geometric fidelity would be likely to 
lose some of the small gradations in image density Information 
sought and their use is not proposed. 

Current and Planned Work 

9. The RBV-1 and -2 images, and the MSS-il have been hand- 
registered and punched for registration pins. A template 
with pins has been made to hold the images during printing, 
and also to locate a series of calibrated 21-step-tablet 
exposures in different locations on the reproduction films, 
in a sequence which will permit more positive control of. 
each phase of reproduction, processing and masking. The 

density step progression will also enable measurements 
to be more closely related to the significant image 
densities, when working with the high net densities noted 
in para 3. 
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10. MasklriK Schematic . Figure 2 is a more complete diagram 
of the various reproduction and masking steps involved 
in the process, and supplants the Schematic shown as 
Figure 2 in the 7 May 1973 Progress Report. 

11. Normalizinp!: image Densities, and Atmospheric Effects . 

The following assumptions have been made in normalizing 
image densities, for the estimation of differences in 
atmospheric effects between the different spectral bands: 

a. Image densities are linear over a 2.0 D range, within 
the 0.4 D to 2.4 D, +.12 D ERTS positive film repro- 
ductions being used for the experiment. 

b. Within the positive. 2.0 D range, 0.0 D represents 100%, 
and 2.0 D 0% sensor response; thus the density of heavy 
white cloud would be 100% response; this in turn is 
checked against the density of the lightest step on the 
15-step tablet printed with the image. The two densities 
should be about equal, deviations being due to local 
exposure and processing differences in the reproduction. 
This was confirmed by measurement. 

c. As noted in (a), the minimum density in the positive 
(100% sensor response) is made high enough - 0.4 D 

+ ,12D - to ensure all image densities are on the 
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FIGURE 2 - Masking Diagram 














straight line of the reproduction curve; if a cloud 
has 0.50 D, this is 100? sensor response, and 0.50 D 
should be subtracted from other image densities to ' 
find their proper relationship to sensor response on 
a 0 % to 100? scale, 

d. However, densities representing less than 20? of full 
scale sensor response are reproduced at gammas less 
than 1,0*. If this effect is plotted with density 
against log^o Percent sensor response, the curve 
shown in Figure 3 results. The dark parts of the 
positive image falling on the curved line have less 
contrast, and less density than they should. In 
effect the curved line shows the non-linearity of 
sensor response. 

The effect of (d) is critical to images of the ocean part 
of the scene, which is generally recorded well below 20? 
sensor response. Because of this, and the difference in 
scene contrast compression caused by the atmosphere between 
the green and red regions, corrective action is necessary 
to adjust image contrast to a level more representative 
of the real values, before images from the two bands can 
be combined for masking; otherwise, spurious effects will 
result. 

* Data Users Handbook, Appendix H 
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% Max 
Sensor 
Output 

Gradient 

(Gamma) 

Normalized 
Y 1.0 Film D • 

Effective 
Film Density 

o 

« 

CM 

0.75 

1.70 

1.27 

3.2 

0.80 

1.50 

1.20 

5.0 

0.85 

1.30 

1.10 

8.0 

0.90 

1.10 

0.99 

12.7 

0.95 

0.95 

0.90 

20 Plus 

1.00 

0.70 

0.70 



Figure 3 - Y 1«0 line represents linear sensor output plotted 

against positive Image density, where 100 ^ response Is 0.0 D 
and 0 ^ Is 2.0 D. Sensor recording gradient (gamma) Is less 
than Y 1*0 between 2 % and 205S, resulting In decreased densities 
and contrast In the darker Image areas, as shown by the curved 
line. Effective density is found by multiplying the theoretical 
Y 1«0 density by the actual local gamma, viz 1.70 D x y 0.75 = 
1,27 D, Average gamma between 4 ^- 10 ^ is 0 . 60 . 
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12, In Figures 4, 5 and 6, Image densities for the same cloud, 

deep water and shallow areas In MSS-4, -5 and -7 have been 
normalized by subtracting the cloud densities as in the 
procedure described in (c) above, the corrected densities 
being plotted on the respective curves, MSS-7 has been 
used as a control for the other bands in evaluating 
atmospheric effects, since these are least in the MSS 
infrared region, and a calm, deep, clear water ocean 
surface should approach as IR absorption will be 

almost complete. In the MSS-7 image, the normalized water 
density should therefore represent the lowest spectral 
reflectance in the scene, and the lowest percentage of 
sensor response, on the assumption that percent scene 
spectral reflectance within any MSS band is equivalent 

to percent sensor response, 

13. The relative percentages of sensor response/scene reflectances 
as viewed and recorded from the sensor through the atmosphere, 
are found by dropping lines from the image subject densities 
to their corresponding abscissa values, as shown in the 
Figures, The following is obtained, percent scene reflectance 
being the antllog of sensor response : 
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MSS 

Band 

i\ 

5 

7 


Log % 

Sensor Response 
Deep Water Shallow 
1.25 tjo 1.32 1.58 to 1.78 

.96 to 1 . 0^1 1.28 to 1.34 

.06 to .30 .36 to .48 


% Scene * 
Reflectance 


Deep Water 
17 . 8 ^ to 21^ 

9 ^ to 11 ^ 

1.255 to 2% 


Shallow 
38 ? to 60^ 

19 ^ to 22^ 

2.355 to 3^ 


* Relative to cloud at 100? 


14. These data show clearly the order of atmospheric effects in 
the three spectral regions. In MSS-7, apparent water 
reflectance is 3? or less in both deep and shallow water. 

In MSS-5 ( 580-680 nm), deep water reflectance apparent to 
the sensor has increased to about 10?, 1% of which is probably 
caused by atmospheric haze, since yellow and red are not 
found in the clear ocean water. Increased reflectivity 
in the shallow water is attributable to red and yellow 
light reflected from the bottom. In the MSS-4 490-600 nm 
spectral region, recorded scene reflectances would be 
predominantly green and yellow. In this area of the 
Caribbean the color of the barren, deep water is blue, 
with a component of green from upwelling subsurface light, 
and reflected skylight not in excess of 3^-5?. The I 8 ? to 
21? apparent ocean scene reflectance is therefore mostly 
contributed by Rayleigh and Mle scattering in the atmosphere. 
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The net haze factor is estimated to be in the order of 
15 /S in MSS-^ and will apply at least equally to the similar 
spectral region recorded by RBV-1, in MSS-5 the haze factor 
is about \s% to 85s, half that of MSS-4, and will apply to the 
RBV-2 image. 


15. Net Effect of Scene Contrast Degradation 

Figure 3 shows an average gamma gradient line drawn on the 
part of the curve between and 10^ sensor response, where 
- it departs from linear y 1»0 recording/reproduction, and 
where the ocean deep water surface will be recorded. Gamma 
in this region is about O.6O, and scene contrast will be 
reduced by that factor accordingly. In Figure 7, 15% and 
7.5% haze curves have been plotted for MSS-4 and -5 images 
respectively, and gamma gradient lines have been drawn on 
them for the to 10% scene reflectance regions. With 
15% haze factor, the real y 1.0 scene contrast is reduced 
to y 0.19, and with 7.5% haze, to y 0.^2. The net effects 
of the gamma reductions are: 


a. In MSS-^ the deep water densities fall on the y 1.0 
linear part of the sensor response curve (Figure 4), 
and only correction for 15% haze is required. The 
final reproduction gamma correcting scene contrast 
in this image would be 


y 1.0 
y .19 


= y 5.26. 
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NORMALIZED DENSITY 



FiKure 7 - Effect of 7.5^ and 15^ estimated haze factors In 
reducing H%-10% scene reflectances from y 1.0 to 
Y 0.42 and Y 0.19 respectively. 
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In MSS -5 the deep water densities fall on the non- 
linear part of the sensor curve, where gamma measures 


(— >• further 

Response 

reduces this to y .^12, the net gamma in this part of 
the image equalling .42 x .75, or y 0.32. In restoring 
scene contrast, final reproduction gamma would then be 


Y 1.0 

Y .32 


, or Y 3.12. 


Provided the density ranges are matched, the two images 
so reproduced would then be equivalent in relative spectral 
scene reflectance and contrast. The density range and 
gamma of the negative image to be used as a mask with any 
positive can then be made to match the positive in opposite 
polarity. 


l6. Error sources . Four main sources of error are inherent 
in this technique; Potential radiometric errors of 10% 
in the RBV and 6% in the MSS; gamma and density variations 
in the input Images, the tolerances being Y + 0.1 and 
locally, 0.06 D to 0.10 D; the accuracy of density measurement 
with the MacBeth Model TD-102 Densitometer, + 0.02 and in 
repeatability + 0.01; and the accuracy of sensitometric 
control in each stage of image reproduction. In addition, 
the careful selection and density measurement of the same 
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part of the image scene in each spectral band must be 
considered. In the set of images being used, comparison 
of the GSPC step tablets indicates that gamma deviations 
do not exceed + 0.05, and the use of the same cloud as 
reference for lOOyS sensor response adjusts the 2.0 useful 

image density range automatically. Many density readings 

i 

on similar image areas in the different bands are consistently 
repeatable to an average of +0.01D. By careful sensitometry , 
the reproduction stages are controlled within the reading 
tolerances of the densitometer, and each stage is plotted 
to confirm that key image densities are reproduced on the 
straight line of the reproduction material. The remaining 
largest unknown variables in the input images are local 
density differences may be 0.06 D to 0.10 D, and their 
radiometric errors. 

17. However, with this particular image set, as shown in 

Figures 4, 5 and 6, the normalized data agree remarkably 
well with theory. It is concluded that in this image set, 
radiometric and local density variations are not large 
enough to drastically Influence results obtained with this 
technique of interactive image processing for data reduction. 
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Ocean Water Color 
Assessment from 
ERTS-1 RBV and MSS 
Imagery, MMC #538 

D. S. Ross PR175 

SIGNIFICANT RESULTS 

1, Interpretation Techniques Development 

/ 

a. Interactive Image Processing ; 

Stationary residual radiometric correction errors In 
RBV-1 and RBV-2 bulk processed 9.5-lnch positive 
transparencies from ERTS-1 were substantially reduced, 
by photographic masking with negative equivalent RBV 
Images taken over a relatively cloud-free ocean area. 
Reported at Symposium on Significant Results Obtained 
from ERTS-1, 5-9 Mar 73 NASA/GSPC (Paper I4) "Masking 
Bulk RBV Images to Reduce Stationary Residual Errors 
In Radiometric Correction" 

b . Interactive Image Processing ; 

An MSS-A negative combined with an RBV-2 positive was 
used to produce a new Image, with the aim of cancelling 
Image densities In a spectral region common to both 
Images. The MSS-4 band records In the 490-600 nm region, 
and the RBV-2, 580 to 680 nm; 580 to 6 OO nm being 
Included In both. This 20 nm part of the spectrum 


21 



appears yellow-orange and is in a region where 
reflectance from sandy bottom in shallow water is 
transmitted. A shallow water area between North 
and South Bimini Islands, where such bottom reflectance 

I 

is strong in the MSS-4 image, is cancelled in the new 
image made by masking. (ERTS-1 E 1007-15165 30 Jul 72) 

c. Effects of the Atmosphere . ERTS-1 bulk processed 

Images are closely held to a specified density range 
and reproduction gamma. By knowing, or estimating 
with reasonable accuracy, the percentage spectral 
reflectance of light and dark subjects in the images, 
a density normalization procedure has been worked out 
for estimating relative effects of the atmosphere in 
degrading scene contrast in the green and red Images 
relative to the Infrared MSS-7. Cloud and calm, 
clear water surfaces are used for light and dark 
references. The solution is largely graphical, and 
can account for changes in scene reproduction gamma at 
sensor response levels less than 20 % of full response. 
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' Ocean Water Color 

Assessment from 
ERTS-1 RBV and MSS 

I Imagery, MMC #538 

j D. S, Ross PR175 

■ SIGNIFICANT RESULTS 

2, Marine Resources & Ocean Surveys 

a. Bathymetry (Alternatively, Interpretation Techniques 

Development, Image Enhancement) 

Photo-optical and electronic density slicing were applied 
to ERTS-1 E 1007-15165-^, 30 Jul 72, an area in the 
Caribbean showing deep ocean water, and shallow areas 
on the Great Bahama Bank ranging from 0.5 meter or less 
to 18 meters. The density slicing processes were adjusted 
to correlate apparent water radiance to bathymetric 
contours shown on C&GS Chart 1112 (1:466,9^0). A number 
of large areas corresponding to water depths of 2 meters 
or less, 5 to 10 meters and 10 to about 20 meters were 
isolated by both processes. Where clear water and 
uniformly reflective bottom was found, clear of marine 
growths, the photo-optical and electronic image density 
slicing processes proved effective in delineating areas 
where depth was in the order of 5 meters, +1 meter. Results 
were reported at Symposium on Significant Results Obtained 
from ERTS-1, 5-9 Mar 73, NASA/GSPC (Paper MI9). 
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